Cavity QED determination of atomic number statistics in optical lattices 
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We study the reflection of two counter-propagating modes of the light held in a ring resonator 
by uhracold atoms either in the Mott insulator state or in the superfluid state of an optical lattice. 
We obtain exact numerical results for a simple two-well model and carry out statistical calculations 
appropriate for the full lattice case. We find that the dynamics of the reflected light strongly depends 
on both the lattice spacing and the state of the matter-wave held. Depending on the lattice spacing, 
the light held is sensitive to various density-density correlation functions of the atoms. The light 
held and the atoms become strongly entangled if the latter are in a superfluid state, in which case 
the photon statistics typically exhibit complicated multimodal structures. 

PACS numbers: 42.50.-p, 42.50.Dv, 42.50.Nn 



I. INTRODUCTION 

The study of ultracold atoms in optical lattices has 
been a very active field of research both experimentally 
and theoretically in recent years. One motivation to 
study these systems is that they provide clean realiza- 
tions of important models of condensed matter physics 
[l| such as the Bose-Hubbard and Fermi-Hubbard mod- 
els 0, Q , spin systems i4!] and the Anderson lattice model 
5| . A particularly well-known example is the prediction 
2| and observation of the Mott-insulator to superfluid 
transition 0| in trapped Rubidium atoms. An im- 
portant advantage of ultracold atoms over conventional 
solid-state systems is that they offer an exquisite de- 
gree of control over system parameters such as interac- 
tion strengths, densities and tunneling rates. Ultracold 
atomic systems are also very versatile in that the atoms 
involved can be either fermionic or bosoniCjand can be 
associated into molecules of either statistics [1, [13, [IH ■ 
Furthermore the dimensionality of the systems can be 
tuned from three to two and one-dimensional. Impor- 
tant other potential applications of ultracold atoms in 
lattices include quantum information [T^ . ITsj and im- 
proved atomic clocks 

The ability to characterize the many-particle state of 
atomic fields is an important ingredient of many of these 
studies. While a full characterization of the field would 
ideally be desirable, a great deal can already be learned 
from the atomic density fluctuations at each lattice site 
and from the intersite density correlations. The count- 
ing statistics of atoms in an optical lattice have previ- 
ously been studied by using spin changing collisions [lH , 
atomic interferences in free expansion 16] and the con- 
version of atoms into molecules via photoassociation [l7j 
or Feshbach resonances [9]. The main goal of this pa- 
per is to propose and analyze an alternative cavity-QED 
based method to measure these properties. 

Our proposed scheme is an extension to trapped ul- 
tracold atoms of the familiar diffraction technique used 
to probe order in crystalline structures. One important 
new aspect is that since we wish to measure the quan- 
tum fluctuations of the atomic density we cannot simply 



use classical radiation, which only provides information 
on some average of the atomic occupation numbers. In- 
stead, we exploit the quantum nature of the light field 
and the fact that it can become entangled with the atoms 
to probe the number statistics of the matter- wave field. It 
is for this reason that a cavity-QED geometry is particu- 
larly attractive: as is well known, high-Q optical cavities 
can significantly isolate the system from its environment, 
thus strongly reducing decoherence and ensuring that the 
light field remains quantum mechanical for the duration 
of the experiment. 

More specifically, we consider two counter-propagating 
field modes in a high-Q ring cavity, coupled via Bragg 
scattering off the atoms trapped in an optical lattice. 
If one of the counter-propagating modes is initially in 
vacuum, this is analogous to a quantum version of the 
Bragg reflection of X-rays by a crystal. 

The remainder of this article is organized as follows: 
Section |TT] describes our model and presents several gen- 
eral results of importance for the following analysis. In 
particular, it shows that the dynamics of the light fleld 
strongly depends on the manybody state of the atomic 
field as well as on the lattice spacing. Section Hill presents 
a series of results for the case of a simple two-well sys- 
tem, analyzing the properties of the Bragg-reflected light 
field for atoms in a Mott insulator state and for a su- 
perfluid described both in terms of a number-conserving 
state and of a mean-fleld coherent state. In particular, 
we find that these two descriptions lead to major differ- 
ences in the properties of the scattered light. Section HVl 
then turns to the case of a large lattice. Finally, section 
|V] is a conclusion and outlook. 



II. 



MODEL 



We consider a sample of bosonic two-level atoms with 
transition frequency uja trapped in the lowest Bloch band 
of a one-dimensional optical lattice with M lattice sites 
and a lattice spacing d, see Fig. [1] We assume for sim- 
plicity that the effects of tunneling and collisions are fully 
accounted for by the initial state of the matter-wave field, 
and neglect them during the subsequent scattering of the 
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Figure 1: (Color online) Atoms trapped in an optical lattice 
with lattice constant d in a ring resonator and interacting 
with two counter-propagating modes with wave vectors ±k. 



weak quantized probe fields, taken to be two counter- 
propagating (plane-wave) modes of a ring cavity with 
mode functions £±k{z), wave- vectors ±fc and frequency 
uJk- We further assume that these fields are far detuned 
from the atomic transition, |A| = |wfc — uja\ 2> j,^Ir, 
with 7 the atomic linewidth and f2/{ the vacuum Rabi 
frequency, so that the excited electronic state of the two- 
level atoms can be adiabatically eliminated. 

Expanding the field operators for the atomic excited 
state ip^'^^z) and ground state ^'^^^z) in the Wannier 
basis of the lowest Bloch band as 



where Cm^ and Cm'' are the bosonic annihilation operators 
for an atom in the excited and ground state at lattice site 
m and Tpm^ and ipm^ the corresponding wave functions, 
this system is described by the effective Hamiltonian 
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(1) 



where we have neglected collisions and tunneling, as al- 
ready discussed, and also ignore all loss processes. In this 
expression 



M-l 

N{d) = J2 e^^'"'^'^c(^)ta(.) 

m=0 



and the coupling frequency is 
2p2 
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(2) 



(3) 



with p the dipole matrix element of the atomic transition. 
Since only the ground electronic state is relevant in this 
model, we drop the superscript {g) in the following for 
notational clarity. 

Our proposed detection scheme relies crucially on the 
observation of the coherent dynamics resulting from the 



Hamiltonian ([T]) , hence it is important that the associated 
characteristic time scales are significantly faster than 
than those associated with losses and with the atomic 
motion in the lattice. To illustrate that this is within 
current experimental reach, we consider explicitly the ex- 
perimental parameters of Klinner et. al. . Assuming 
that the atoms are confined to much less than an opti- 
cal wavelength at the antinodes of the lattice potential, 
we can approximate the coupling constant between the 
cavity modes and a single atom as 
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where V is the mode volume of the cavity and p w 2.3 x 
lQ-29 jg |-j^g relevant dipole matrix element for the 
D2 line of ^^Rb. For a detuning of A = —1 GHz from 
that transition, a cavity length of 0.1 m, and a mode 
waist of 100 /im we find g « 50 s^^. 

As we show later on, the fiuctuations of the opera- 
tor N{d) are central in the determination of the spectral 
properties of G. As these typically scale like the square 
root of the total number of atoms we find that for a sam- 
ple of 10^ atoms the relevant frequencies for the coherent 
evolution are of the order of 50 kHz. Such frequencies 
are several times larger than the decay rate of state-of- 
the-art optical cavities with large mode volume, 17 kHz 
in the experiments of Klinneris et. al. The spontaneous 
emission from the atomic excited state is negligible at 
this detuning. Indeed, the fact that in the experiment 
of Klinner a splitting of the normal modes has been ob- 
served is direct experimental evidence that the coherent 
dynamics can be made dominant over the losses. 

Turning now to the characteristic time scales for the 
center-of-mass atomic dynamics, specifically the interwell 
tunneling rate and the two-body collision rate, they can 
be controlled with high accuracy, respectively through 
the depth of the lattice potential and via a magnetic 
Feshbach resonance. It is therefore possible to operate 
under experimental conditions such that both are much 
smaller than the coupling strength between atoms and 
light. This justifies "freezing" the atomic dynamics re- 
sulting from tunneling and collisions once the probe fields 
are turned on. 

It is well known that, depending upon the ratio be- 
tween tunneling and two-body collisions, bosonic atoms 
in an optical lattice can undergo a transition from a su- 
perfiuid to a Mott insulator state. With this in mind, we 
consider initial atomic states that correspond to these 
two extreme situations, specifically: (a) a Mott insulator 
state with a well-defined atom number in each well, (b) a 
state where each well is in a coherent state, a mean-field 
approximation of the superfluid state, and (c) a more re- 
alistic description of atoms in the superfluid state with 
a flxed total number of atoms N. These three states are 
given explicitly by 



IV'Mott) = \no, ni,. . . , UM-i), 



SFl) 



|ao,ai, 



(4) 
(5) 
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where |amp is the mean number of atoms in well m, and 




where M — %/ NIM^ is a normalization constant. 

We further assume for simplicity that the mode prop- 
agating in the ~k direction along the ring is initially in a 



vacuum while the other mode is in a Fock state with ntot 
photons. Our results are qualitatively independent of the 
exact state of the light field, but we will point out the 
modifications brought about by a coherent state instead 
of a number state when appropriate. 

The intensities of the two counter-propagating modes, 
(a^Ofc) and (alj,a_fc), are easily obtained from the solu- 
tion of the Heisenberg equations of motion for a±k(t), 



a-k{t) 
ak{t) 



cosG{d)t 
-igt((f) sin G{d)t 



-iQ{d) smG{d)t 
cos G{d)t 



a-fc(O) 
afc(O) 



(7) 



where 



III. 



DOUBLE- WELL POTENTIAL 



N{d) 



Q{d) = - 

^N{d)N{-d) 

is unitary and we have introduced the operator 



(8) 



G{d) = g\lN{d)N{-d) 



(9) 



describing the coupling between the two modes. The 
closed form solution ([7]) follows from the observation that 
in the absence of tunneling the optical field only couples 
to conserved quantities of the atomic field. Note however 
that the number statistics of the atoms at each site is 
conserved for many cases of practical interest even in the 
presence of tunneling. This is for instance the case for 
the ground state of the Bose-Hubbard model. Our results 
are valid for such cases as well. 

Expanding the atomic states in the number states ba- 
sis jrioj ^ii • ■ ■ J nM-i), where the operators G and Q are 
diagonal, we then find for the intensity of the Bragg- 
reflected field 



ntot 



E 



no,ni,...,nAf-l 

X COS 2a;(no, ni, . . . , nM-i)i, 



(10) 



where P„ 



...^riM-i is ttis probability to find atoms 
in the m'^ well and LoijiQ^ni, . . . ,nM-i) is the corre- 
sponding eigenvalue of G. 

As we see in the following sections, this expression ac- 
counts implicitly for the well-known property that Bragg 
scattering depends strongly on the well separation d. In 
particular, the light scattered from neighboring wells in- 
terferes either constructively or destructively for well sep- 
arations equal to A/2 or A/4, respectively, two situations 
relatively easy to realize experimentally. 



Before considering the situation of a general optical lat- 
tice, we discuss the simpler case of a double- well poten- 
tial, where explicit analytical results are readily obtained. 
Despite its simplicity, this case already exhibits many of 
the properties of the full M-site lattice and hence pro- 
vides valuable intuition for its understanding. We first 
discuss the Bragg-rcflected intensity for well separations 
d — X/2 and d = A/4, and then turn to an arbitrary well 
separation and to the analysis of the photon statistics. 



A. Reflected intensity 

1. Well separation d — X/2 



From Eq. (|T0|) , most properties of the reflected inten- 
sity can be understood from an analysis of the eigenval- 
ues of the operator G{d) describing the coupling between 
the forward- and backward-propagating light fields. For 
d A/2 we find from the definition (O of N{d) that 

G(A/2) - .g(no + ni), 

indicating that the intensity of the Bragg-scattered light 
is only sensitive to the total number of atoms N , a signa- 
ture of the constructive interference of the fields reflected 
off the two wells. Specifically, 



[n-k 



ntot 



Pn sin^ gNt 

N=0 



(11) 



where P/v is the probability that the total atom number 
in the two wells is N. 

Both the Fock state |'0Mott) and superfluid state |'0sf2) 
have a well defined total number of atoms so that the light 
field simply undergoes harmonic oscillations between the 
+k and —k directions. These two situations are illus- 
trated in the upper and lower parts of Fig. [21 
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Figure 2: Reflected intensities (n_fe(t)) for well separation 
d = \/2 and for atoms (a) in a Fock state |'i/>Mott), (b) in a 
superfluid state in mean field approximation Ii/jsfi) and (c) 
in the number-conserving superfluid state |?/'sf2). The mean 
number of atoms in each well is nine in all three cases. 



In contrast, if the atomic field is described by coherent 
states for each well, the total number of atoms is uncer- 
tain and consequently the reflected intensity is comprised 
of oscillations at the eigenfrequencies associated to all 
possible combinations of atom numbers, leading to col- 
lapses and revivals similar to those familiar from the two- 
photon Jaynes-Cummings model [H,[23|, see the middle 
curve in Fig. [51 The Bragg-reflected intensity periodi- 
cally collapses and revives after characteristic times 
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and 



^revival — 1 9 1 



(12) 



(13) 



where A(no + "-i) — \/ (("-o + "-i)^) — ('"■o + ni)'^ is the 
standard deviation of the total number of atoms. Equa- 
tion (I12p is important in that it indicates that the inten- 
sity of Bragg-scattered light is indeed sensitive to fluc- 
tuations in the atom number. As such, it offers a clear 
signature to distinguish the mean-field description from 
the number-conserving description of the atomic super- 
fluid state. 



2. Well separation d = X/ 4, 

The light reflected from the two wells is now tt out 
of phase and hence interferes destructively. Mathemati- 
cally, this is reflected by the coupling between the two 
counter-propagating modes being now proportional to 
the difference between the atom numbers in the two wells, 

G(A/4)=g|no-ni|, 
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Figure 3: Reflected intensities {n_fc(f)) for well separation 
d = A/4 and for the atoms (a) in Fock states with equal atom 
numbers, (b) in coherent states of equal mean atom number 
and (c) in the number-conserving superfluid state. The mean 
number of atoms in each well is nine in all cases. In (c) the 
solid line is for A'' = 18 and the dashed line is for A'^ = 17. 



SO that the intensity of the backscattered light fleld is 
given by 



{n-k{t)) 
ntot 



P\AN\sin^ g\^N\t, (14) 

|A7V|=0 



where aat] is the probability that the atomic population 
difference between the two wells is AiV or — AiV. 

If the atomic fields in both wells are in Fock states, the 
sum in Eq. (fH)) reduces to just one term and the light 
intensity undergoes sinusoidal oscillations at frequency 
g\nQ — ni\. If the atom numbers are equal the two optical 
modes become decoupled, see Fig. Ela). 

If each well is in a coherent state, on the other hand, 
the number difference is not known with certainty and 
accordingly the oscillations of the intensities of the two 
modes contain frequency components corresponding to 
all possible number differences. In the case of equal 
mean atomic numbers in the two wells the reflected in- 
tensity rises to half the initial intensity of the forward- 
propagating fleld, remains constant at this value and 
drops to zero at time Trevivai- The dips in the reflected 
intensity can be thought of as 'remnants' of the collapse 
and revival features exhibited by the d = X/2 case. Sim- 
ilarly to that case the half-width of the dips is deter- 
mined by the uncertainty in the atom number differ- 
ence, which is half of the uncertainty in the total number 
of atoms '22'], so that the initial rise time is given by 

^collapse = 2T^^ji^p^^^ = 1/5V(«0) + (m). 

Finally, for atoms in the number-conserving superfluid 
state |'0SF2), odd differences in atomic numbers are im- 
possible for an even total number of atoms, hence all 
possible oscillation frequencies differ by at least 2g (or ig 
for the oscillations of the field intensities). For an even 
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Figure 4: Construction of the eigenvalues of G for d — A/10 
and probability of occurrence of each frequency for the two 
wells in coherent states with coherent amplitudes ai = 02 — 
3. 

total number of atoms the revival time is cut in half, 
see Fig. [3jc). For odd numbers of atoms the reflected 
intensity is equal to ritot at times t = Tr/2g, 37r/2(7, . . . 
Apart from these times the reflected intensity behaves 
very similar to the coherent state case. The collapse time 

collapse collapse ' 

3. Arbitrary well separation 

The expressions (|lip and for the backscattered 
photon numbers in the two cases considered so far indi- 
cate that the marginal probability distributions Pn and 
-P|AAr| can be directly determined from the reflected in- 
tensity by means of Fourier transforms. Unfortunately, 
the complete number statistics Pno,ni for uq atoms in 
the first well and ni atoms in the second well cannot be 
reconstructed from these marginals. 

The dynamics of the two optical modes becomes more 
involved for arbitrary separations, as the eigenvalues of 
the operator G are no longer integer multiples of g but 
in general irrational numbers that depend on the occu- 
pation numbers no and fii in a nontrivial manner. One 
consequence is that the perfect recurrences observed for 
d = X/2 and d = A/4 no longer occur and the dynamics 
of the backscattered intensity is typically non-periodic. 

The backscattered intensity is determined by the spec- 
tral properties of the operator G as we have seen, and 
these can be understood by means of a simple geomet- 
ric construction: For a number state with uq atoms in 
the first well and ni atoms in the second well the eigen- 
value of G is given by the hypothenuse of the triangle 
obtained by drawing a straight line of length gng and 
then continuing with another straight line of length gni 
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Figure 5: Backscattered intensity for atoms in state IV'sfi) 
with nine atoms in each well on average and a well separation 
d = A/10. 



at an angle 2kd with respect to the first line. This con- 
struction is shown in Fig. 01 which also illustrates the 
probability of each frequency occurring if each well is in 
a coherent state with coherent amplitude uq = ai = 3. 
The dots along the tlq line are drawn in shades of gray 
with darker shades being more likely and lighter shades 
being less likely. The ni lines originating from each dot 
are drawn with the same convention, and the probability 
for a given combination of no and ni, obtained by mul- 
tiplying the probabilities for no and ni, is again shown 
in gray scale. Hence the most likely eigenfrequencies and 
the spread in frequencies can be read off as the distances 
of the final dots from the origin. A typical example of 
Bragg-refiected intensity is shown in Fig. [5l This exam- 
ple is for nine atoms per well in the coherent state IV'SFi) 
and a well separation d = A/ 10. 

Figure [6] shows the spectrum of G as a function of d 
for state |'0Sfi) with = a\ = 3. Each frequency is 
again weighted with its likelihood, and the eigenfrequen- 
cies have been collected in bins of the size of the pixels in 
the figure. The special cases of d — A/ 2 and c? = A/4 are 
easily recognized, as for those well separations only inte- 
ger multiples of g can occur. For general well separations, 
the spectra lack any obvious structure. 

The number-conserving superfluid state |'0SF2) corre- 
sponds roughly to a "diagonal cut" through the con- 
struction in Fig. 2] such that ng decreases by one if n\ 
is increased by one. Accordingly there are now much 
fewer possible frequencies as can be seen in Fig. [T] 
Furthermore, the width of the frequency distribution is 
no longer independent of d. We have A|no — ni| = 

y^f (l-cos2fcd) [13. 



B. Photon number statistics 

A major advantage of the two-well case over the full 
lattice case is that we can easily solve the Schrodinger 
equation for the coupled atom-cavity system numerically. 
This integration is greatly facilitated by the conservation 
of the number of atoms in each well and the total number 
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Figure 6: Spectrum of G as a function of d. The darkness 
of each eigenfrequency corresponds to its likelihood. In this 
example both wells are in a coherent state with an average of 
9 atoms per well. 
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Figure 7: Spectrum of G as a function of d for atoms in state 
|i/'SF2) with 18 atoms. 



of photons, which we exploit by integrating in the invari- 
ant subspaces of constant no, hi and fik + h-k- From the 
full solution we can gain more detailed information about 
the dynamics of the photons than the average reflected 
intensity considered so far. 

Here we consider the number statistics of the Bragg- 
reflected photons, 



Pn^, = (|n_fe)(n_fe|). 



(15) 



The calculation of Pn_f, involves a trace over the well oc- 



cupation numbers. As a consequence the number statis- 
tics is, in complete analogy to the reflected intensity in 

Eq. pUj) . the sum of the number statistics Pn"^'"^'* for 
each set of occupation numbers (rto,"i), each weighted 
with its probability of occurring in the atomic state, 



p — \^ p p(no,' 



Figure 8: Photon number statistics Pn_k{t) for each well in 
a Fock state with no = 9 and ni = 10 atoms, a well separa- 
tion of d = A/4 and 10 photons initially in the +k direction. 
Lighter shades of grey correspond to higher probabilities. 
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Figure 9: (Color online) Photon number statistics for atoms 
in state \1psF2) for A*' = 18 atoms and a well separation d — 
A/4. Also indicated are the elementary sinusoidal oscillations 
corresponding to number differences one and two. 



Figure [5] shows one of the building blocks Pi"°'"'^(t) 
for well occupation numbers uq — 10 and ni — 9 with 
d — A/4. In this case the two counter-propagating modes 
are linearly coupled to each other, with coupling fre- 
quency g\nQ — ni\. Figure [5] shows the number statistics 
for atoms in state |'0sf2)- This example is representative 
of the typical situation, with non-trivial photon statistics 
and 'complicated' dynamics of the light field that cannot 
be adequately described by means of the reflected in- 
tensity only. As the elementary probability wavepackets 
Pnl"k"^\t) dephase with respect to each other the pho- 
ton statistics evolves through complex patterns, with a 
large uncertainty of the photon number. The reflected 
light intensity exhibits large fluctuations the description 
of which requires all moments of the intensity up to order 

«tot- 

There are however special times when several of the 
elementary photon statistics building blocks add up to 
yield peaks at specific photon numbers. Prominent such 
times are t — 7T/2g^^ and t — ir/Ag^^ for d ~ A/4. 
At these times all odd frequency components (i.e. odd 
number differences) lead to — ntot while all even 
number differences lead to h-k = 0. This is illustrated 
in Fig. [TO] for the atomic state iV'SFi)- For the fairly 
large mean number of atoms per well considered in that 
example the probabilities for even and odd frequencies 
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Figure 10: (Color online) Photon number statistics at times 
(a) t — -R /Ag~^ and (b) t — -n /2g~^ for atoms in state IV'sfi) 
with average atom numbers no = n\ = Q for well separation 
d = A/4. 



are approximately equal so that the probabilities to find 
all photons reflected or no reflection at all are also nearly 
equal. For atoms in state |'0SF2) the photon statistics 
at time t — TT/Ag~^ look very similar to Fig. [TUTb). 
It is important to emphasize that the light field does 
not evolve into a 'NOON' state (|ntot,0) + |0, ntot))/^^, 
which would show the same behavior. Rather, the light 
field is entangled with the atoms in the Schrodinger cat 
state 



m = TT/2g)) = iVsc) 

^|ntot,0) [|ao,(-ir°'ai) + |-ao, 



-(-l)"-ai)] + 10, ntot) [\ao, (-l)"-ai) 



|-ao,-(-ir-ai>] 



which is reflected by the atomic Q-function in the first 
well, 

Q(a) = (^sc|a>(a|^sc) = ^ (e-l"-"«l' + e-l"+"°l') . 

For d = X/2 the number statistics looks qualitatively 
similar to the A/4 case. The most likely frequencies are 
however typically higher because of the constructive in- 
terference of the light reflected off the two wells that leads 
to the coupling frequency being g{ho + fii). 

In practice, Fock states are very hard to realize. Typ- 
ically the +k mode will start out in a superposition of 
states with different photon numbers such as a coherent 
state. As we have mentioned, this case is qualitatively 
similar to the Fock state case and here we would like to 
explain in more detail why this is so. Because the sys- 
tem Hamiltonian conserves the total number of photons 
fi-k + 'T-fc, no interference between states with different 
total photon numbers can arise in any observable that 
preserves the total photon number as well. The reflected 
intensity and the photon number statistics are clearly 
observables of this type. For a general state the expecta- 
tion value of these observables is the sum of the expec- 
tation values for the various photon numbers weighted 
with the probability of that photon number. The conse- 
quence is simply that one has additional fluctuations in 
the reflected intensity due to the uncertainty of the pho- 
ton number in the initial state. Another way to see that 
the intensity behaves the same for all initial states of the 
+k mode is to insert the formal solution for the modes 
of the light field Eq. ([7]) into the expectation value of 
the intensity. The terms containing operators of the —k 
mode vanish and we end up with the expectation value 
of the intensity in the +k mode at time t = multiplying 
the trace of the atomic operators over the atomic states. 



IV. OPTICAL LATTICE 
A. Well separation d = A/2 

This section extends the results of the double well anal- 
ysis to the case of an optical lattice with a large number 
of wells. The case of a well separation c? = A/2 is again 
particularly simple since in that case G is simply propor- 
tional to the total number of atoms in the lattice, 

G = gY,n,n = gN. (16) 

m 

Both the Fock state and the number-conserving super- 
fluid state \'4'SF2) have a well-defined total number of 
atoms N, hence the intensity of the optical field under- 
goes simple sinusoidal oscillations at frequency gN be- 
tween the two counter-propagating modes. 

The situation is more complicated in the mean-field de- 
scription of the superfluid state, in which case the matter- 
wave field at each lattice site is in a coherent state with 
mean number of atoms rim = {N)/M. In the limit of a 
large lattice, M ^ 1, the central limit theorem permits 
to approximate the probability distribution of the total 
number of atoms to an excellent degree as 

Pn = , exp . (17) 

[ 2(Ar) J 

In the limit of large {N) this distribution is sharply 
peaked at (iV) and the dynamics of the light is largely 
characterized by oscillations at frequency g{N) between 
the two counter-propagating modes. However the finite 
variance of the total atomic number distribution results 
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in a collapse of these oscillations on a time scale 

1 



rpX/2,SFl,latt. 
collapse 




(18) 



The photons undergo roughly 2y (N) oscillations before 
this collapse. Note that since the allowed frequencies are 
discrete and integer multiples of g, the oscillations also 
show complete revivals, with a revival time given by Eq. 



B. Well separation d — X/A 

The situation is markedly different for a well separation 
d = A/4. In that case G is proportional to the difference 
in the total number of atoms trapped on even and odd 
lattice sites, 



G^g 



J2 f'rn- J2 



m even 



odd 



g\Ne - No 



(19) 



and the light field probes fluctuations in the occupation 
differences of neighboring sites. This is a useful property 
since fluctuations of this type, while being typically dif- 
ficult to measure, can serve to distinguish various many- 
body quantum states such as the Mott insulator state 
and the superfluid state. Note that if all lattice sites 
are in Fock states with equal atomic population, then G 
vanishes identically for an even number of lattice sites 
and is very small when the number of lattice sites is odd. 
In typical experiments the atoms are also subject to a 
harmonic trap and the added trapping potential leads to 
shells with constant occupation number separated from 
each other by edges with a superfluid component. Equa- 
tion (|19p shows that the reflected light can serve as a 
good detector of these edges. 

For the atoms in the mean- field superfluid state IV'SFi) 
the coupling strength can again be analyzed using the 
central limit theorem. To this end we use the fact that 
Ne and —No are the sums of a large number of indepen- 
dent random variables with means {N)/2 and —{N)/2, 



respectively, and standard deviation y {N)/2. From the 
resulting Gaussian probability distributions for N,, and 
No we readily find the probability distribution for N^—No 
by using that the sum of two Gaussian random variables 
is again a Gaussian. For the intensity of the reflected 
light it is actually |iVe — No\ that is important, see Eq. 
(fT4|) . The probability distribution for this quantity is 
readily found from those for N^^o as 



-P|7Ve-Af„| = (2 - Sn^^nJ 



1 



: exp 



'2t:{N) 

where d is the Kronecker delta. 



{Ne - No)' 
2{N) 



(20) 




Figure 11: (Color online) Probability distribution of the cou- 
pling frequencies (eigenvalues of G for d = A/4 with atoms 
in state \1psF2) for 10 lattice sites with 20 atoms. The bars 
are the exact probabilities determined numerically and the 
diamonds show the distribution obtained from Eq. (|2ip . The 
line is an aid to the eye. 



The analysis of the state \1psF2) requires a bit more 
attention. Since Ng is the sum of many operators it is 
natural to assume that its distribution is Gaussian, and 
hence fully characterized by its mean (V'SF2|-^e|'0SF2) = 
N/2 and standard deviation AA^e — VN /2. Since every 
choice of Ng automatically determines No = N — N^ the 
probability distribution for |A^e — -^o| is 



V2 



: exp 



[Ne - No) 



21 



2N 



(21) 



Figure [TT] shows an example of the numerically ob- 
tained exact probability distribution and compares it to 
the approximate result of Eq. (|2ip for as few as 10 lat- 
tice sites and 20 atoms. Clearly, the two results already 
agree quite well even for these modest numbers, so that 
the expression (I^Tj) can be considered as essentially ex- 
act for the much larger numbers of lattice sites and atoms 
typically encountered in practice. 

Note that if one assumes that the occupation numbers 
of the lattice sites contributing to N^ are independent 
when calculating the probability distribution for N^ with 
the central limit theorem, one finds that the standard 
deviation of the resulting distribution is too large by a 
factor of \/2. In this sense, fluctuations of the number 
difference between even and odd sites are suppressed in 
the state |V'sf2)- 



C. General well separation d 

The calculation of the possible eigenvalues of G and 
their probabilities of occurring for the atomic states 
I'f/'SFi) and |7/'sF2) by means of the central limit theorem 
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can readily be extended to well separations d = q\/2p, 
where q and p are integers. The lattice sites can be 
grouped into p distinct classes according to the roots of 
unity e^'^'^ = e^^'^'/P with l = OA,...,p~l. The total 
number of atoms in each class |2JI , 



M/p-l 

= ^ ni+mp 
m— 

enters the frequency with the common coefficient e^'^'^^lv 
and each class comprises many lattice sites, so that the 
central limit theorem can be used to calculate the asso- 
ciated number distribution P^q^ 25]. The frequencies are 
then obtained as a function of the few macroscopic oc- 
cupation numbers iVj, which are independent to a good 
approximation. Thus the probabilities for the frequencies 
uj{Nq, 7Vi, . . . , 7Vj,_i) = g\ YXZl iVze2*"'/f I are simply the 
product of the probabilities for each macroscopic occupa- 
tion Pjv, . These probabilities can however not be written 
down in closed form, since the frequencies are no longer 
sums of Gaussian variables. 

If p is large or d is an irrational fraction of A, the 
behavior of the system becomes identical for atoms in 
states IV'SFi) and |'0sf2) in the limit of a large number 
of lattice sites. In this case we can think of iV((i) as 
the result of a random walk of M steps of average stride 
length lyN) /M in the complex plane and the frequency is 
proportional to the distance from the origin to the final 
point. The distribution of possible frequencies ui becomes 
quasi-continuous with 



2uj 



exp 



(22) 



independently of whether the atoms are in state IV'SFi) or 
|V'SF2)- This is illustrated in Fig. [12] for atoms in state 
IV'SFi) and in Fig. [T3] for atoms in state \ipsF2)- It is 
remarkable that the numerically obtained true frequency 
distribution is in a good agreement with the approximate 
Eq. (I22p even for the low number of atoms and lattice 
sites considered here, and even more so, that this agree- 
ment is already obtained for p — 5. As one might expect, 
we find much better agreement between the actual fre- 
quency distributions and Eq. (|22p for the atom and site 
numbers of Figs. [T^ and [T51 if the lattice spacings corre- 
spond either to a larger p or to irrational fractions of the 
optical wavelength, as is the case in the inset of Fig. [T^ 
which is for d = -\/2A/10. 

One can calculate the reflected intensity in closed form 
from the frequency distribution (1^^ by approximating 
the sum in Eq. (jlOp by an integral, 



{n-k) 
ntot 



/•OO 

/ diLjP{uj) sin^ uit 
Jo 



2VM 



te- 



erfi 



9{N)t 



M 



(23) 




Figure 12; (Color online) Probability distribution of frequen- 
cies for atoms in the mean-field superfiuid state IV'SFi) with 
ten lattice sites and ten atoms for d = A/10. The solid line is 
the distribution Eq. (|22|l . and the histogram the result of an 
exact numerical diagonalization of G{d). 



3, 
Q_ 




Figure 13: (Color online) Probability distribution of fre- 
quencies for atoms in the number-conserving superfiuid state 
|'!/'SF2) with ten lattice sites and 18 atoms for d = A/10. The 
solid line is the distribution Eq. (|22p and the histogram a re- 
sult of the numerical diagonalization of G{d). The inset shows 
the same distribution for d — 



where erfi(a;) 



JJ^ dze~^ is the complex error func- 
tion. The intensity converges to (ri_fe) = ntot/2 after a 
transient of duration TcoUapse ~ "^/^/(si-^))- 



V. CONCLUSION 

In this paper we have investigated the Bragg reflection 
of a quantized light field off ultracold atoms trapped in 
an optical lattice as a function of the lattice spacing. 
We have considered atoms in a Mott state as well as a 
superfiuid state, the latter described both in the mean- 
field approximation and in a number-conserving form. 
We have studied both a simplified two-mode model that 
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lattice 


atomic 


J.-1 
collapse 


spacing 


state 


2 wells 


lattice 




h/'Mott) 


U 


U 


7 A 


h/'SFl) 


2ff\/{n-o) + (ni) 


2^7^/w 




h/'SF2) 










I'i/'Mott) 








d - A 


l^/'SFl) 


g\J{no) + (Jii) 


3\/l-2/7ry^{iV) 




h/'SF2) 




g^/l - 2/nVN 




IV'Mott) 








general d 


IV'SFl) 


2ff\/{n-o) + (n-i) 


Vl-^/4ff(7V)/VM 




\lpSF2) 


2gy^f (1 - cos2/cd) 


^/l - TT/4gN/VM 



Table I: Collapse times for the various lattice spacings and 
atomic states. See the comments in footnote [21] regarding 
the validity of the double well results for l-i/isFi) and \ipsF2)- 



allows us to carry out most of the calculations explicitly 
or numerically as well as the full lattice case, where our 
analysis mostly relies on statistical arguments. 

Our results show that the dynamics of the light 
field strongly depends on the many-particle state of the 
matter-wave field and thus can serve as a diagnostic tool 
for that state. Depending on the lattice spacing the light 
field can probe fluctuations of various density correla- 
tion functions on the lattice. These fluctuations typi- 
cally lead to collapses of the oscillations of the reflected 
intensity. From the standpoint of probing the atomic 
many-particle state the lattice spacing d = A/4 seems the 
most promising, since in that case the three atomic states 
give rise to dramatically different reflection signals: The 
Mott state leads to perfect sinusoidal oscillations, the su- 
perfluid state described in the mean-field approximation 
to collapses with revivals after t = ir/g and the number- 
conserving superfluid state to revivals after t — i: /2g~^ 
for even total atom numbers and anti-revivals at the same 
times for odd total atom numbers. Also, the sensitivity 
of the reflected intensity to the number differences on 
neighboring sites means that it can be used as a detector 



for the edges of the density plateaus in the Mott insulator 
state in a harmonic trapping potential. 

The double-well case can easily be handled by numeri- 
cally solving for the coupled dynamics of atoms and light 
field. The solution shows that the atoms and light field 
evolve into non-trivial entangled states if the atoms are 
in a superfluid state. For those states the photon number 
statistics is typically many-modal and the characteriza- 
tion of the state of the light field merely in terms of the 
reflected intensity is incomplete. 

Most results obtained for the double-well case imme- 
diately carry over to the case of a lattice, where we 
have shown that for general well separations, the system 
behaves the same for the mean-field and the number- 
conserving superfluid state. 

Future work will extend this model to include cavity 
losses thereby making it possible to model more realisti- 
cally the measurement of the reflected photons. It will 
be interesting to determine into which atomic states the 
matter- wave field is projected as a result of this mea- 
surement. Another open question is the determination 
of the best possible choice of lattice period d, or even of 
an optimal sequence of lattice spacings for the purpose of 
reconstructing the full counting statistics of the atomic 
field. We plan to investigate this point in more detail by 
studying the information content of the reflection signal 
by means of a Bayseian analysis. 
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